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THE  LOWER-HYBRID-DRIFT  INSTABILITY  IN  NON-ANTIPARALLEL  REVERSED 

FIELD  PLASMAS 


I.  INTRODUCTION 

An  important  set  of  problems  in  plasma  physics  consists  of 
understanding  the  physical  processes  which  can  occur  in  reversed 
field  plasmas.  For  example,  how  can  magnetic  field  energy  be 
rapidly  converted  into  particle  energy?  How  does  the  topology 
of  the  magnetic  field  change?  Under  what  conditions  can 
topological  changes  occur?  An  enormous  amount  of  work  has  been 
devoted  to  these  (and  related)  questions  over  the  past  20  years. 
Generally  speaking,  this  research  has  been  focused  on  the  inves¬ 
tigation  of  magnetic  field  reconnection  processes.  Some  parti¬ 
cular  topics  of  interest  have  been  field  line  annihilation  (e.g., 
ID  Sweet-Parker  models);  forced  reconnection  (e.g.,  2D  Petschek 
model) ;  and  tearing  instabilities.  In  an  astrophysical  context, 
this  research  has  been  relevant  to  a  variety  of  space  phenomena 
such  as  solar  flares,  interplanetary  D  sheets  and  geomagnetic 
substorms. 

One  process  which  can  be  important  in  reconnection  physics 
is  plasma  microturbulence.  In  general,  plasma  microinstabilities 
can  often  produce  anomalous  transport  of  particles,  momentum  and 
energy.  This  can  be  critical  to  reconnection  process,  especially 
in  the  null  region  (B  *•  0) ,  since  it  allows  the  plasma  to 
"decouple"  from  the  magnetic  field.  Moreover,  anomalous  trans¬ 
port  effects  can  greatly  enhance  the  rate  of  energy  conversion 
from  the  magnetic  field  to  the  plasma.  A  variety  of  micro¬ 
instabilities  which  can  lead  to  fine  scaled  turbulence  have  been 

analyzed  to  determine  their  relevance  to  reconnection  (see 
Manuscript  submitted  July  17, 1981. 


Papadopoulos  (1979)  for  a  review) .  However,  a  drawback  of  many 
of  these  analyses  is  the  assumption  of  a  one-dimensional  magnetic 

A 

field  (i.e.,  B  =  B  (x)e  );  an  assumption  not  usually  justified  in 
space  plasmas.  Other  components  of  the  magnetic  field  leads  to 
the  following  effects.  First,  a  component  of  B  normal  to  the 

A  A 

neutral  line  (i.e.,  B  =  B..e„  +  B„e_)  introduces  field  line  cur- 
vature.  This  effect,  which  is  generally  incorporated  in  micro¬ 
instability  theories  via  an  artificial  gravity,  leads  to  addition¬ 
al  particle  drifts  and  can  be  either  a  stabilizing  or  destabiliz¬ 
ing  influence  depending  upon  the  plasma  conditions.  Secondly, 
a  component  of  B  parallel  to  the  current  (i.e.,  B  =  ^e^  +  Bzez) 
introduces  magnetic  shear.  Physically,  this  corresponds  to  the 
situation  where  the  magnetic  fields  on  either  side  of  the  neutral 
are  not  anti-parallel.  Magnetic  shear  leads  to  Landau  resonances 
of  particles  and  waves,  and  the  coupling  of  cross-field  modes  to 
parallel  propagating  modes.  Magnetic  shear  is  generally  a  stabi¬ 
lizing  influence  on  instabilities.  Since  the  magnetic  shear  in¬ 
duced  Landau  resonances  and  parallel  mode  couplings  are  strongly 
dependent  on  spatial  position,  the  analysis  of  this  effect  re¬ 
quires  a  nonlocal  theory.  In  contrast,  the  magnetic  curvature 
induced  particle  drifts  are  not  strongly  dependent  on  spatial 
position  and  can  be  accurately  analyzed  with  a  local  theory. 

In  this  paper  we  analyze  the  effect  of  magnetic  shear  on  the 
lower- hybrid-dri ft  instability  and  discuss  some  of  the  implica¬ 
tions  of  magnetic  shear  as  it  regards  the  dynamics  of  reversed 
field  plasmas  in  the  magnetosphere.  He  choose  the  lower-hybrid- 


drift  instability  since  it  is  the  roost  likely  instability  to  be 
excited  in  reversed  field  space  plasmas  of  interest  (i.e. ,  the 
earth's  magnetopause  and  magnetotail  (Huba  et  al.,  1978;  Gary 
and  Eastman,  1979)).  Moreover,  the  anomalous  transport  properties 
associated  with  this  instability  can  be  important  to  reconnection 
processes  (Huba  et  al. ,  1977;  Drake  et  al. ,  1981;  Huba  et  al., 
198Q).  We  focus  here  on  the  effect  of  magnetic  shear  on  the 
lower- hybrid-drift  instability  since  it  has  been  shown  to  have 
a  strong  stabilizing  effect  on  this  mode  (  Krall,  1978].  We  will 
not  consider  here  the  effect  of  field  line  curvature  since  it  has 
a  much  weaker  influence  on  the  lower-hybrid-drift  instability 
(Krall  and  McBride,  1977;  Rajal  and  Gary,  1981). 

A  self-consistent  theory  of  the  lower-hybrid-drift  instabi¬ 
lity  in  a  sheared,  reversed  field  plasma  is  difficult  to  develop. 
Major  complications  arise  because  of  finite  plasma  8  and  electron 
temperature  effects  (i.e.,  electromagnetic  coupling,  electron  VB 
drift-wave  resonance,  finite  electron  Larmor  radius  effects). 

The  inclusion  of  these  effects,  within  the  context  of  a  weakly 
nonlocal  stability  analysis,  leads  to  three  coupled,  partial 
differential  equations  involving  complex  velocity  integrations. 

At  this  early  stage  of  analysis,  we  do  not  attempt  such  a  com¬ 
prehensive  theory.  Rather,  we  use  a  simpler  theory  here  to  ex¬ 
amine  the  dominant  effect  of  magnetic  shear  on  the  lower-hybrid- 
drift  instability  with  emphasis  on  understanding  the  implications 
of  this  mode  for  reversed  field  plasma  dynamics.  Specifically, 
we  simplify  the  analysis  by  using  electrostatic  theory  with 
cold  electrons.  However,  we  do  include  an  Appendix  which 


discusses  the  effect  of  electromagnetic  coupling  and  electron  VB 
drift-wave  resonances  in  the  shear  stabilization  criterion.  A 
more  comprehensive  theory  will  be  presented  in  a  future  report. 

The  scheme  of  the  paper  is  as  follows.  In  Section  II  we 
present  a  general  discussion  of  the  lower-hybrid-drift  instabi¬ 
lity  and  the  effect  of  magnetic  shear  on  it;  showing  specific 
results  of  our  analysis,  i.e.,  a  marginal  stability  curve  show¬ 
ing  the  stabilizing  influence  of  shear.  We  also  apply  these  re¬ 
sults  to  a  particular  sheared,  reversed  magnetic  field  configu¬ 
ration.  In  Section  III  we  discuss  the  implications  of  these  re¬ 
sults  on  reconnection  processes  in  the  earth's  magnetosphere. 
Finally,  an  analysis  incorporating  finite  8  effects  is  presented 
in  the  Appendix. 
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D.  THEORY 

A.  Assumptions  and  Plasma  Configuration 
The  plasma  configuration  and  slab  geometry  we  consider  is 
described  as  follows.  The  ambient  magnetic  field  is  in  the  y-z 

A  A 

plane  (B  =  Byey  +  Bzez)  an<*  is  only  a  function  of  x.  The  density 
is  also  a  function  of  x  but  we  take  the  ion  temperature  to  be 
finite  and  constant.  For  simplicity,  we  take  the  electron  tem¬ 
perature  to  be  zero  (i.e.,  T  /T.  -*-0).  Finite  electron  temperature 
effects  are  discussed  in  the  Appendix.  Equilibrium  force  bal¬ 
ance  on  an  ion  fluid  element  in  the  x  direction  requires  viy~v^i 

where  =  (v^2/2£K)3  6m  n/3x  is  the  ion  diamagnetic  drift  velo- 

k 

city.  Here,  v^  =  (2T^/m^)  is  the  ion  thermal  velocity  and 
=  eBQ/m^c  is  the  ion  Larmor  frequency.  We  can  relate  the 
ion  diamagnetic  velocity  to  the  mean  ion  Larmor  radius  and  scale 
length  of  the  density  gradient  by  =  rLj/2Ln  where 

=  vi/fii  and  Ln  =  ^  n/3x)  .  The  electrons  are  assumed  to 

be  magnetized,  while  the  ions  are  treated  as  unmagnetized.  This 
is  reasonable  since,  in  treating  the  lower- hybrid-drift  instabi¬ 
lity,  we  are  considering  waves  such  that  <<  w  <<  and 
k2r^2  >>  1.  Only  electrostatic  oscillations  are  considered 
(electromagnetic  coupling  is  discussed  in  the  Appendix)  and  we 
assume  that  the  plasma  is  weakly  inhomogeneous  in  the  sense  that 
rLe2  Sm.  n/3x)2  <<  1  and  rLe2  (3  6m  B/3x) 2  <<  1.  We  assume  that 
the  local  equilibrium  magnetic  field  is 

B(x)  *  B0(xQ)(e2  +  (x-x0)/Lg4y) 
in  the  vicinity  of  x  =  x  (i.e.,  (x-x  )/L  <<  1)  where 

sJ  O  D 


(i) 


L  =  (3$/3x)  *and  4>  =  tan  *(£  B_).  Thus,  L_  is  the  scale  length 
s  y  z  s 

characterizing  the  magnetic  shear. 

If,  for  the  moment,  we  consider  B  =  BQez,  then  the  plasma 

configuration  just  described  is  unstable  to  the  kinetic  lower- 

.  k 

hybrid-drift  instability  when  1  >VH./v.  > (m  /m.)  (Davidson  et 

Cll  X  6  X 

al. ,  1977).  The  instability  is  driven  by  the  cross-field  cur¬ 
rent  and  is  excited  via  an  ion-wave  resonance  (i.e.,  inverse 
Landau  damping) .  The  waves  are  characterized  at  maximum  growth 
by  a)r  -v  kyVdi  *  u>£h  ,  y  ^  wr  ,  ky  pes  -  1  and-k •  B  =  0  where 
Pes  =  (T^/me)  For  modes  such  that  k*B  ^  0  (i.e.,  k|.  ^  0), 

electron  Landau  damping  reduces  their  growth  rate  or  stabilizes 
them,  depending  on  the  magnitude  of  k^  .  We  now  let  B  =  By(x)ey  + 
Be  as  in  Eq.  (1) ,  which  introduces  magnetic  shear  as  shown  in 
Fig.  1.  The  magnetic  field  rotates  in  the  y-z  plane  as  a  func¬ 
tion  of  x  so  that  kii  is  also  a  function  of  x.  At  x  =  x  we  note 

II  o 

that  k jj  =  0  (k*B  =  0)  but  at  x  =  x^,  k|j  f  0  ( jk * B  ^  0).  Thus,  the 

dispersive  properties  of  the  plasma  are  also  a  function  of  x 

(e.g.,  electron  Landau  damping  and  parallel  mode  coupling  can 

occur  at  x^  but  not  at  xq) .  Making  use  of  Eq.  (1),  we  use  the 

prescription  kz (x)  =  kZQ  +  k^  (x-xq)/Ls  to  incorporate  magnetic 

shear  into  the  analysis  where  we  choose  k  =  k  (x=x  )  =0. 

zo  z  o 


B.  Dispersion  Equation 

Within  the  context  of  the  assumptions  outlined  in  the  pre¬ 
vious  sub-section,  the  equation  which  describes  the  lower-hybrid- 
drift  instability  in  a  sheared  magnetic  field  is  given  by 
[Davidson  et  al.,  1978]. 


6 


(2) 


where 


a2* 

A  +  [B-C(x-x  )]S  =  0 


dx 


k  A  =  - 


*V  L‘ 


B  = 


—el  r 

k2Vi2  L 


k2p  2-  z"  U,) 

es  i 


kV  - . 


] 


Pi  ,2  2  di  .  . 

*—!>  k  p„ - - —  +  1  + 


eS  U> 


(3) 


(4) 


C  =  J2*-  -L- 

u,2  L  2 
s 

and  “'po2  =  4line2/n>a,  v±2  =  2Ti/mi,  Peg2  =  Hj/»e°e2>  * 

=  (v^2/2^i)d  Bm.  n/dx  and  =  (ut-k^V^^) /kv^ .  In  writing 

Eq.  (2)  we  have  assumed  k2X  2  <<  1  and  ui  2  »  £J  2  where 

D  pe  e 

,  2  _  2 /0  2 
XD  "  vi  /2u,pi  * 

Equation  (2)  is  in  the  form  of  Weber's  equation  and  the 
eigenf requency  is  defined  by 

B  =  (2m  +  1)  (AC)3* 


(S) 


(6) 


where  m  is  the  mode  number  (i.e.,  m  =  0f  1,  2,  ...).  The  branch  is 
chosen  according  the  outgoing  energy  prescription  of  Pearlstein 
and  Berk  (1969).  The  associated  eigenfunction  is 

<J>  =  <t>0  Hm(ux)  exp  [-a2x2/2]  (7) 

where  a  -  (C/A)^  and  Hm  is  the  Hermite  polynominal  of  order  m. 

In  general,  a  numerical  analysis  is  required  to  solve  Eq.  (6) 
(using  Eqs.  (3)-(5)J.  However,  we  first  consider  a  limiting  case 
to  obtain  an  analytical  solution  to  Eq.  (6) . 
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r  * 


We  consider  the  weak  drift  limit,  Vdi  «  vif  so  that  the  Z 
functions  in  Eqs.  (3)  and  (4)  can  he  approximated  by  Z  ( )  “*  i/rr. 
The  dispersion  equation  then  becomes 


BC.W  -  1  ♦  ^2Pes2  -  ^ 


+  i 


“pe  _1 _ 

“  kyLs 


(2m  +  1) 


(8) 


The  first  imaginary  term  is  a  destabilizing  term  due  to  inverse 
ion  Landau  damping.  The  second  imaginary  term  is  the  stabilizing 

effect  of  magnetic  shear.  It's  origin  in  the  analysis  is  a  term 

kii 

*  (~w  2  t*ie  ma9neti-zed  electron  response.  Physically,  magnetic 

shear  leads  to  stabilization  since  it  allows  wave  energy  to  propa- 
cfate  away  from  the  excitation  region  (i.e.,  where  kjj  =  0). 

The  real  frequency  is  given  by 


•r  “  kvdi''tl  + 


(9) 


where  shear  corrections  to  have  been  neglected.  The  mode 
is  stabilized  when  Im  D(iu,k)  =  0  or 


„  .2  2 

/?  Vdi  k  pes _ 

(2m+1)  vi  (l+k2p  2)  2 

6  S 


(10) 


where  Eq.  (9)  has  been  used  and  the  subscript  cr  refers  to  the 

critical  value  of  L  /L  .  The  maximum  value  of  the  RHS  of  Eq. 

n  s 

2  2 

(10)  occurs  for  k  p  =  1.0  so  that  all  wavenumbers  are  stable 

\  r  gg 


when 


( 2m+l ) 


for  > 


/? 

8~ 


(2m+l) 


(ID 
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Note  that  the  higher  order  inodes  (m  jt  0)  are  more  easily 

stabilized  by  shear  than  the  lowest  order  mode  (m  =  0)  (i.e., 

L  (m  jt  0)  >  L  (m  =  0)), 
s  s 

We  now  relax  the  weak  drift  assumption  and  solve  Eq.  (6) 

numerically.  Figure  (2)  is  a  plot  of  (L  /L  )  vs.  V../V,,  for 

n  s  tn  i 

2  2 

u)pe  /^e  =  100,  0  =  0.0,  m  =  0  and  T0  =  0 .  The  growth  rate  is 

maximized  with  respect  to  kpeg  so  that  the  maximum  value  of 

(L  /L  )  necessary  to  stabilize  the  mode  is  obtained.  Two 
n  s  cr 

regions  are  labeled:  stable  (L  /L_  >  (L _/L_)  )  and  unstable 

ns  ns 

(L  /L  <  (L  /L  )  ).  For  V,./v.  =  1.0  we  find  that  (L  /L) 

ns  ns  cii  i  ns  qj- 

0.325  which,  surprisingly,  is  in  good  agreement  with  our 
analytical  result  [Eq.  (11)]  even  though  the  drift  is  not  weak. 

C.  Application  to  Reversed  Field  Plasmas 

In  order  to  apply  these  results  to  a  reversed  field  plasma, 
we  consider  the  following  magnetic  field  profile 

B  =  B  [sin  e„  +  cos  *■  tanh  £  e  l  (12) 

o  z  y  z  kz 

Equation  (12)  describes  a  reversed  field  plasma  (i.e.,  the  z- 
component  reverses  direction  at  x  =  0)  with  the  field  undergoing 
a  total  directional  change  tt  +  0  in  the  y-z  plane.  Thus,  if 
0=0  then  the  standard  one-dimensional  Harris  profile  is 
recovered  (i.e.,  anti-parallel  field)  with  a  discontinuous 
rotation  at  x  =  0.  On  the  other  hand,  for  0^0  the  field  con¬ 
tinuously  rotates  in  the  y-z  plane  and  remains  finite  every¬ 
where.  This  profile  (Eq.  (12) )  is  a  rough  approximation  to  the 
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Fig.  2  —  Plot  of  (L^/Lglcr  vs.  V^/Vj  for  Wpe2/ne2  =  100  and  Te  =  0.  The 
lower-hybrid-drift  instability  is  stable  (unstable)  for  L-/Lg  >  (L^/Lglg  (L„ /L* 
<  (Ln  /Lg)^).  The  curve  is  dashed  for  V^/Vj  <  (me/mi)'^  *  0.16  because  the  ions 
are  magnetized  in  this  regime  and  our  theory,  strictly  speaking,  is  not  applicable. 


magnetic  field  of  the  nose  of  the  magnetosphere  of  the  earth. 

It  also  approximates  the  distant  geomagnetic  tail  when  no  normal 
field  component  is  present. 

We  plot  Ln/Ls  vs.  x/A  for  several  values  of  0  (0  =  5°,  10°, 

15°,  20°,  25°,  30°,  45°)  in  Figure  3.  We  note  that  L /L  varies 

li  s 

2-3  orders  of  magnitude  over  the  range  of  x/A  shown.  Also,  as  6 

increases  so  does  Ln/Lg,  as  anticipated.  In  the  region 

x/a  ^  0.25,  the  shear  is  strong  (i.e.,  L  ?>  L  )  and  the  theory 

S  it 

discussed  in  this  paper  is  not  adequate.  Thus,  we  restrict 
our  attention  to  regions  such  that  Ln  •>  Lg . 

We  combine  the  results  of  Figs.  (2)  and  (3)  in  Fig.  4 
which  plots  o  vs.  x^/X  for  several  values  of  Wvt  <vdi/vi  -  °-25 
0.50,  1.00  which  correspond  to  A/rT  .  =4.0,  2.0,  1.0,  respectively)- 
Here,  Xp/A  represents  the  linear  penetration  distance  of  the 
lower-hybrid-drift  instability.  That  is,  we  expect  the  modes  to 
be  stable  for  x  <  Xp  because  of  shear  stabilization.  We  find 
that  (1)  as  the  ion  diamagnetic  drift  increases  (i.e.,  the  cur¬ 
rent  sheet  becomes  thinner)  the  mode  can  penetrate  closer  to  the 
"null"  region  (i.e.,  x  =“  0)  and  (2)  as  0  increases,  which 
increases  the  magnetic  shear,  the  penetration  distance  Xp/A 
becomes  larger.  Nonlocal  analysis  of  the  lower-hybrid-drift 
instability  in  a  field  reversed  plasma  (Huba  et  al.,  1980)  has 
found  that  the  dominant  mode  is  localized  at  a  position 
x/a  “•  1.25  for  Tg  =  0.  Thus,  even  for  0  =  45°  the  fastest 
growing  mode  is  not  stabilized  due  to  shear.  However,  higher 
order  modes,  which  penetrate  closer  to  x  =  0,  are  expected  to  be 


0.25  0.50  0.75  1.00  1.25  1.50 

Xp/X 


Fig.  4  -  Plot  of  0  vs.  XpA  for  Vdi/vj  =  1.00,  0.50,  0.25.  Here,  0  represents  the 
strength  of  the  magnetic  shear  in  a  reversed  field  plasma  (i.e.,  0*0:  no  shear, 
fields  anti-parallel;  0  i  0:  shear,  fields  non-antiparallel)  and  xpA  is  the  linear 
penetration  distance  of  the  lower-hybrid-drift  instability  to  the  neutral  line.  The 
mode  is  linearly  stable  (unstable)  in  the  region  x  <xp  (x  >xp).  As  the  amount  of 
shear  increases,  i.e.,  0  becomes  larger,  the  mode  is  stabilized  further  away  from 
the  neutral  line  (x  =  0). 
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affected  by  shear.  These  results  are,  however,  highly  profile 
dependent.  If  the  dominant  shearing  of  the  field  occurred  in 
the  region  where  the  LHD  was  localized,  rather  than  in  the  null 
region  (the  situation  implied  by  Eq.  (12)),  then  the  stabilizing 
effect  of  shear  would  be  much  more  pronounced. 


01.  DISCUSSION 

The  purpose  of  this  brief  report  has  been  to  study  the  in¬ 
fluence  of  magnetic  shear  on  the  lower-hybrid-drift  instability 
in  reversed  field  plasmas.  As  mentioned  in  the  introduction, 
microturbulence  associated  with  such  modes  can  play  an  important 
role  in  the  dynamic  evolution  of  a  reversed  field  plasma  via 
its  anomalous  transport  properties.  The  lower- hybrid-dri ft  in- 
stabiltiy,  although  linearly  stable  in  the  field  reversal  region 
near  the  neutral  line  (Huba  et  al.,  1980),  does  play  a  dramatic 
role  in  the  evolution  of  an  anti-parallel  field  reversed  plasma 
(Drake  et  al.,  1981;  Winske,  1981;  Tanaka  and  Sato,  1981).  It 
has  been  shown  both  theoretically  and  by  computer  simulations 
that  the  mode  causes  magnetic  flux  to  diffuse  towards  the  neutral 
line  which  leads  to  an  enhanced  current  density  at  the  neutral 
line.  Eventually  microturbulence  penetrates  to  the  null  region 
and  permits  field  line  reconnection/annihilation  to  occur,  there¬ 
by  dissipating  magnetic  energy.  In  the  case  where  the  reversed 
field  is  non-antiparallel  (i.e. ,  the  field  is  sheared)  we  have 
found  here  that  the  lower-hybrid-drift  instability  will  be 
linearly  unstable  further  away  from  the  neutral  line  than  the 
case  where  the  field  is  completely  anti-parallel  (Fig.  4).  Thus, 
magnetic  shear  has  an  inhibiting  effect  on  the  penetration  of 
the  lower-hybrid-drift  mode  toward  the  neutral  line.  We  anti¬ 
cipate  that  this  means  that  the  evolution  of  a  non-antiparallel 
reversed  field  plasma  will  differ  from  that  of  an  anti-parallel 
one.  For  example,  if  the  shear  is  sufficiently  strong  the  mode 


may  take  substantially  longer  to  penetrate  to  the  neutral  line 
or  may  not  penetrate  at  all.  This  conjecture,  of  course,  is 
based  on  our  linear  analysis  and  must  be  substantiated  by  further 
(nonlinear)  analysis.  Also,  we  reemphasize  that  the  results 
presented  in  this  paper  are  not  comprehensive  since  we  have 
neglected  important  finite  8  and  Te  effects  (i.e.,  coupling  of 
electromagnetic  and  electrostatic  fluctuations  and  VB  resonances). 
Nonetheless,  our  results  are  qualitatively  correct  although  they 
can  be  improved  quantitatively.  We  are  presently  developing 
a  more  comprehensive  theory  for  finite  8  and  T@  plasmas  and 
will  report  our  results  in  a  future  publication. 

Two  regions  of  the  magnetosphere  in  which  the  magnetic 
field  can  reverse  direction  (in  one  component)  and  is  also  sheared 
are  the  nose  and  magnetotail.  In  both  these  regions  reconnection 
processes  are  believed  to  occur  and  can  be  important  to  the  dyna¬ 
mic  interaction  of  the  solar  wind  and  the  magnetosphere.  In  the 
case  of  the  nose,  the  angle  between  the  incoming  IMF  and  the 
earth's  geomagnetic  field  varies  from  0  to  it,  A  thin  magneto¬ 
pause  boundary  layer  exists  (whose  width  is  roughly  rT . )  over 
which  the  B  field  undergoes  a  directional  change.  In  the  magne¬ 
totail  it  is  known  that  a  substantial  crosstail  magnetic  field 
(By  *  15y)  can  exist  at  times  [Akasofu  et  al.,  1978] .  Such  a 
magnetic  field  can  introduce  a  strong  shear,  i.e.,  the  magnetic 
field  undergoes  a  rotation  as  one  passes  from  the  north  to  south 
lobe. 

Crooker  (1978)  has  suggested  that  the  nose  reconnection  only 


occurs  in  regions  where  the  INF  and  geomagnetic  field  are  anti¬ 
parallel  (no  shear) .  Classical  reconnection  theories  predict 
reconnection  even  if  the  fields  are  non-antiparallel  [ Ugai, 

1981]  so  that  seme  "anomalous"  process  may  be  responsible  for 
inhibiting  reconnection  in  this  instance.  Moreover,  for  recon¬ 
nection  events  observed  at  the  nose,  the  time  scale  for  the 
energization  of  the  plasma  appears  to  be  slower  for  more  strongly 
sheared  cases  (non-antiparallel  fields)  than  for  non-sheared 
fields  (anti-parallel  fields)  (Russell,  private  communication) . 

We  suggest  one  possibility  for  such  an  effect  is  the  influence 
of  magnetic  shear  on  the  lower- hybrid-drift  instability  (or  other 
instabilities)  and  its  associated  microturbulence. 
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APPENDIX 


! 

j 


We  derive  an  expression  for  (L  /L  )  which  includes  finite 

n  3  cr 

3  effects  (i.e.,  electrostatic-electromagnetic  coupling,  elec¬ 
tron  drift-wave  resonance) .  In  order  to  make  the  analysis 

tractable,  we  consider  the  limit  3,  «  1  and  T  <<  T. .  The  non- 

i  e  i 

local  dispersion  equation  is  given  by  [Davidson  et  al.,  1977] 


2 

A  +  [B-C(x-x  )  2  ]  <f) 
dx2  ° 


0 


(Al) 


where 


A 


(A2) 


B 


C 


(A3) 


(A4) 


where  u>pa2  =  47mea2/ma,  va2  =  2Ta/ma,  peg2  =  (^/n^)/^2. 

Pi  =  PunT^/B2 ,  Vdi  =  (vi2/2Q±)b£n  n/Sx,  sfi  =  v.2/ve2,  PA  *  (■- 

kvd^)/kv^,  vf2  ■  2Qeour/kyefi  and  efi  =  &  Hn  B/3x.  The  finite  3 

corrections  included  in  (Al)  are  as  follows.  First,  the  coupling 

of  electromagnetic  and  electrostatic  perturbations  arises  from 

2  2  2 

terms  proportional  to  u)pe  /c  k  .  The  3  dependence  can  be  seen  by 


I 


1 


2  2  2  2  2 

noting  that  k  pgg  (“>pe  /c  k  )  -  &^/2.  Second,  non-resonant 
electron  7b  effects  are  contained  in  the  term  proportional  to 
@i  (kV’di/u')  •  This  is  a  fluid-like  response  of  the  electrons  due 
to  the  inhomogeneous  magnetic  field.  Finally,  resonant  electron 
Vb  effects  are  contained  in  the  final  term  of  B  in  Eq.  (A2) . 
Equations  (A2)-(A4)  are  similar  to  those  derived  by  Davidson  et 
al.  (1977)  but,  in  addition,  contain  electron  resonance  terms 
(Huba  and  Wu,  1977) . 

We  now  assume  ^vr/^e  »  1  which  corresponds  to 

0.  <.<  k2p  2  ~  1  and  rewrite  Eqs.  (A2)  -  (A4)  as  follows 
l  es 


(A6) 

(A7) 


The  dispersion  equation  is  again  given  by  Eq.  (1) .  The  real 
frequency  is  the  same  as  Eq.  (9)  but  with  peg  replacing  peg 
where  p  =  p_ _/(l  +  p,/2)  .  To  lowest  order  in  0.,  the  criti- 

6S  6S  X  1 

cal  shear  length  is 


20 


1 


l+k2p  2 
es 


seexp(-se) 

2 A  2 
1+k  "es 


(A8) 

where  sg  can  be  written  as 


s 

e 


2_ 


l+k2p 


es 


(A9) 


Maximizing  (L  /L  )  with  respect  to  wavenumber  yields 
n  S 


✓tt 

4 


exp(-l/Pe) 

(28i)J* 


(AiO) 


where  kp _ =“  kp _ =  1. 

es  es 

Two  interesting  points  concerning  Eqs.  (A7)  and  (A9)  are 
the  following.  First,  the  finite  8^  dependence  in  the  first 
term  of  Eqs.  (A7)  and  (A9)  arise  from  the  electromagnetic  cor¬ 
rection  due  to  6A||  (i.e.,  the  transverse  magnetic  field  fluctua¬ 
tions)  .  The  influence  of  this  correction  is  to  increase  the 
amount  of  shear  necessary  to  stabilize  the  mode  (Davidson  et 
al.,  1978).  That  is,  as  3-  increases  then  the  shear  length  L_ 
necessary  for  stabilization  decreases  so  that  the  mode  is 
harder  to  stabilize.  Physically  this  occurs  because  the  fluc¬ 
tuating  electric  field  associated  with  6A|j  inhibits  free  stream¬ 
ing  electron  flow  along  the  magnetic  field  which  reduces  the 
rate  at  which  energy  can  be  convected  away  from  the  localization 
region  (Pearlstein  and  Berk,  1969) ,  Secondly,  the  second  term 
in  Eq.  (A8)  represents  the  resonant  Vb  correction  which  is  a 
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damping  effect.  This  term  tends  to  decrease  the  amount  of 
shear  necessary  to  stabilize  the  mode.  Thus,  the  finite  B  cor¬ 
rections  have  different  influences  on  the  shear  stabilization 
criterion.  Loosely  speaking,  electromagnetic  effects  are  de¬ 
stabilizing  (i.e.,  a  stronger  shear  is  needed  to  stabilize  the 
mode  from  the  0=0  situation)  while  the  resonant  VB  effects 
are  stabilizing  (i.e.,  a  weaker  shear  is  needed  to  stabilize 
the  mode  from  the  0=0  situation) .  As  to  which  effect  dominates, 
a  more  careful  analysis  is  needed  which  we  are  presently  develop¬ 
ing  and  will  be  presented  in  a  future  report. 
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